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ABSTRACT 

We analyze the exact general relativistic exact integro-differential equation 
of radiative transfer describing the interaction of low energy photons with a 
Maxwellian distribution of hot electrons in gravitational field of a Schwarzschild 
black. We prove that due to Comptonization an initial arbitrary spectrum of 
low energy photons unavoidably results in spectra characterized by an extended 
power-law feature. We examine the spectral index by using both analytical 
and numerical methods for a variety of physical parameters as such the plasma 
temperature and the mass accretion rate. The presence of the event horizon 
as well as the behaviour of the null geodesies in its vicinity largely determine 
the dependence of the spectral index on the flow parameters. We come to 
the conclusion that the bulk motion of a converging flow is more efficient in 
upscattering photons than thermal Comptonization provided that the electron 
temperature in the flow is of order of a few keV or less. In this case, the 
spectrum observed at infinity consists of a soft component produced by those 
input photons that escape after a few scatterings without any significant 
energy change and of hard component (described by a power law) produced 
by the photons that underwent significant upscattering. The luminosity of the 
power-law component is relatively small compared to that of the soft component. 
For accretion into black hole the spectral energy index of the power-law is 
always higher than one for plasma temperature of order of a few keV. This 
result suggests that the bulk motion Comptonization might be responsible for 
the power-law spectra seen in the black-hole X-ray sources. 

Subject headings: accretion — black hole physics — radiation mechanisms: 
nonthermal — Compton and inverse Compton — radiative transfer — relativity 
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— stars: neutron — X-rays: general 
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1. INTRODUCTION 

Do black holes interact with an accretion flow in such a way so a distinct observational 
signature entirely different from those associated with any other compact object exists? In 
other words can the existence of a black hole be solely inferred from the radiation observed 
at infinity ? 

These are the crucial questions where theoreticians and observers are confronting 
nowadays. Even though there have been now accumulated the enormous observational 
evidences in favor of the existence of black holes still it is fair to say that their existence has 
not been established. Perhaps the proof their existence would have been a much easier task 
if, for instance, an argument would have been advanced, which would: 

a) single out the generic component (or components) of a black hole which is responsible 
for shaping up the unique observed feature associated with black holes 

b) prove that indeed this generic component always results in the same observable 
feature independent of the environmental conditions that the black hole finds itself. 

The luck of such an argument may be traced in the plethora of various accretion 
fiows: accretion in a state of free fall, optically thin or optically thick, accretion disks 
with or without relativistic corrections, shocked flows etc. Of course, this diversity of 
accretion models is highly justified. On physical grounds one expects accretion fiows 
describing a solar-mass black hole accreting interstellar medium to be distinct from those 
fiows describing accretion onto a black hole in a close binary system or from a supermassive 
black hole at the center of an AGN. Viewed from this angle, the detectability of a black 
hole appears to be a rather frustrating issue since it is not clear a priori what type of the 
existing accretion models (if any) would describe a realistic accreting black hole. 

In the present paper we shall show that may be not the case. The distinct feature of 
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black hole spacetime, as opposed to the spacetimes due to other compact objects is the 
presence of the event horizon. Near the horizon the strong gravitational field is expected 
to dominate the pressure forces and thus to drive the accreting material into a free fall. 
In contrast, for other compact objects the pressure forces are becoming dominant as their 
surface is approached, and thus free fall state is absent. We argue that this difference is 
rather crucial, resulting in an observational signature of a black hole. Roughly, the origin 
of this signature is due to the inverse Comptonization of low energy photons from fast 
moving electrons. The presence of the low-energy photon component is expected to be 
generic due, for instance, to the disk structure near a black hole or to Bremstrahlung of the 
electron component from the corresponding proton one. The boosted photon component is 
characterized by a power law spectrum, and is entirely independent of the initial spectrum 
of the low-energy photons. The spectral index of the boosted photons is determined by 
the mass accretion rate and the bulk motion plasma temperature only. A key ingredient 
in proving our claim is the employment of the exact relativistic transfer describing the 
Compton scattering of the low-energy radiation field of the Maxwellian distribution of fast 
moving electrons. 

We will prove that the power law is always present as a part of the black hole spectrum 
in a wide energy range. We investigate the particular case of a non-rotating Schwartzchild 
black hole powering the accretion, leaving the case of a rotating black hole for the future 
analysis. 

The presence of the power law part in the upcomptonized spectra was rigorously 
proven by Titarchuk & Lyubarskij 1995 (hereafter TL95). There it has been demonstrated 
that for the wide class of the electron distributions the power law is the solution of the full 
kinetic equation. 

The importance of Compton upscattering of low-frequency photons in an optically 
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thick, converging flow has been understood for a long time. Blandford and Payne were the 
first to address this problem in a series of papers (Blandford & Payne 1981 and Payne & 
Blandford 1981). 

In the first paper they derived the Fokker-Planck radiative transfer equation which 
took into account photon diffusion in space and energy, while in the second paper they 
solved the Fokker-Planck radiative transfer equation in the case of the steady state, 
spherically symmetric, super-critical accretion into a central black hole with the assumption 
of a power-law flow velocity v{r) (X and neglecting thermal Comptonization. For the 
inner boundary condition they assumed adiabatic compression of photons as r ^ 0. Thus, 
their flow extended from r = to infinity. They showed that all emergent spectra have a 
high-energy, power-law tail with index a = 3/(2 — /?) (for free fall /3 = 1/2 and a — 2), 
which is independent of the low-frequency source distribution. 

Titarchuk, Mastichiadis & Kylafis (1996) paper (hereafter TMK96 and see also the 
extended version in TMK97) presents the exact numerical and approximate analytical 
solutions of the problem of spectral formation in a converging flow, taking into account the 
inner boundary condition, the dynamical effects of the free fall, and the thermal motion of 
the electrons. 

The inner boundary has been taken at finite radius with the spherical surface being 
considered as a fully absorptive. 

Titarchuk, Mastichiadis & Kylafis (1996) have used a variant of the Fokker-Plank 
formalism where the inner boundary mimics a black-hole horizon; no rclativistic effects 
(special or general) are taken into account. Thus their results are instructively useful 
but they are not directly comparable with the observations. By using the numerical and 
analytical techniques they demonstrated that the extended power laws are present in the 
resulting spectra in addition to the blackbody like emission at lower energies. 
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Zane, Turrola, Nobilli & Erna (1996) presented a characteristic method and the 
code for solving the radiative transfer equation in differentially moving media in a curved 
spacetime. Some applications, concerning hot and cold accretion onto nonrotating black 
holes were discussed there. 

In our paper the full relativistic treatment is worked out in terms of the relativistic 
Boltzmann kinetic equation without recourse to the Fokker-Planck approximation in either 
configuration and energy space. 

The relativistic transport theory was developed by Lindquist (1966). He presents the 
appropriate radiative transfer in the curve spacetime. For completeness of our paper we 
delineate some important points of that theory related with the application to the radiative 
transfer in the electron atmosphere. 

We demonstrate that the power-law spectra are produced when low-frequency photons 
are scattered in the Thomson regime (i.e. when the dimensionless photon energy z^E' /meC^ 
measured in the electron rest frame satisfies z' <^ 1). 

The eigenfunction method for the Comptonization problem employed in this paper has 
been offered and developed by Titarchuk & Lyubarskij 1995 (hereafter TL95). 

Giesler & Kirk 1997 have extended the TL95 treatment by accommodating an arbitrary 
anisotropy of the source function. Their results for the spectral index confirm those of TL95 
over a wide range of electron temperature and optical depth; the largest difference they 
found is 10%, occurring at low optical depth. 

The spectral indices related with the eigenvalues of the problem are determined as 
functions of the optical depth of the accreting matter. Thus, for the first time we are able 
to solve the full Comptonization problem in the presence of the bulk and thermal motions 
of electrons. 
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In § 2 and Appendix A we will give the details of the derivation of the general 
relativistic radiative kinetic equation. Section 3 ( and some details in the Appendix B) 
presents the method. We describe the method of separation of variables, and the reduction 
of the whole problem to the specific eigen-problem in the configuration space. We propose 
the numerical solution of this problem by using the iteration method (e.g. Sunyaev 
& Titarchuk 1985) with integrating over characteristics (the photon trajectories in the 
presence of Schwarzschild black hole background). Finally, we summarize our work and 
draw conclusions in §4. 

2. THE MAIN EQUATION 

We begin with considering, background geometry, described by the following line 
element: 

ds^ = - fdf + y + r'dil^ (1) 

where, for the Schwarzschild black hole, / = 1 — r^/r, = 2GM/c^, and t, r, 6, (f are the 
event coordinates with dfl^ = dO'^ + sin^ 9d(p^. G is the gravitational constant and M is the 
mass of a black hole. 

In order to describe the photon radiation field we shall employ the concept of the 
distribution function N. The distribution function N{x, p) describes the number dN 
of photons (photon world lines) which cross a certain spacelike volume element dV at 
x{t,r,6,ip), and whose 4— momenta p lie within a corresponding 3-surfacc clement dP in 
momentum space. It is desirable to choose dV and dP to be coordinate-invariants. Thus 
dN would be invariant as well and the same would be true of N(x, p). 

In Appendix A we present the detailed derivation of the relativistic radiative transfer 
equation expressed through the distribution function N{x, p) and the interaction density 
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function S{N) (see the definition of this function after Eq. A13). 

We will describe the electron component by a local Maxwellian distribution (e.g. 
Landau & Lifshitz 1980, and Pathria 1970 ) 

F{r,Pe)dP = ^-^^''^^^dP (2) 

where is the normalization constant. 

One has to interpret F{r, Pe){—Pgna)dPdV similar to the one implied by (A8), with 
the sole exception that considerations are restricted on the electron phase space. For our 
purpose an arbitrary electron momentum state Pg , can be represented in the form: 

Pe={ J ^ , I '^'"^ I, (3) 

— * 

where the "the thermal three velocity" V stands for a convenient parameterization of the 
electron phase space. We will take j3 = meC^/kTe, while ti^ stands for the hydrodynamical 
four velocity of the inflowing plasma which may be represented relative to the local 
orthonormal frame in the form 



with the negative sign in takes into account the convergent nature of the fluid flow. 

Note that, as a result of the hydrodynamic bulk motion, the local Maxwellian distribution 
exhibits a coupling of the thermal velocity V with the hydrodynamic bulk motion v, and 
one gets 



1 + cos6'— 



C2 



(5) 



We will discuss the coupling effect in §4 and we will consider this issue in detail in our next 
publication. 
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Within a 4— volume dW at the event x there is a decrease in the original number of 
world lines due to absorption and scattering out of momentum range dP, given by (cf. the 
right hand side of Eq. A 13) 

-k{x, p)n{x)N{x, p) dW dP. (6) 

Here n{x) is the proper number density of the electrons interacting with the photons, 
namely the number density of electrons as measured in their own local rest frame, and 
k{x, p) is the invariant absorption coefficient or invariant opacity. The k— opacity is related 
to the usual scattering cross-section Ug via expression (see Lindquist 1966) 

K = E-Gs. (7) 

On the other hand, there is increase due to pure scattering out of all other 4— momentum 
ranges dP' into dP given by 

n{x) dW dP j dP'K(x, p)C(x; p' ^ p)N(x, p'). (8) 

Thus the transition probabihty ({x;p') can be expressed in terms of the differential cross 
section dag/ {dEdO) (see Lindquist 1966) as 

«(.,p')Cfep'-p)^f^. (9) 

Taking into account only Compton scattering of photons off the background electrons, one 
may covariantly write the transfer equation (see Eq. A13) in the following form: 

DN r 

= J Nir, P')k{x, p')ax; p' ^ p)dP' 

-N{r, P) j N{r, P')k{x, p')ax; p ^ p')dP'. (10) 

The first term in the right hand side describes the increase in the photon world lines over 
the infinitesimal phase space cell centered around P while the second term describes the 
processes of depletion. 
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Recall that the scattering cross-section of photon from an electron in the electron's rest 
frame is described by the Klein-Nishina formula 



-UeCTT 
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(11) 



:i+ a [1+^(1-0] J 

where z — hv/me(? is a dimensionless photon energy, ^ is the cosine of scattering angle, 
gt is the Thomson cross— section and that all quantities in the right hand of the above 
formula are computed in the rest frame of the electron one may explicitly write the transfer 
equation on the black hole background. 

By rewriting (10) for the orthonormal frame of (1), [see Eq. (A25)] we get the following 
equation: 
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(12) 



The scattering kernel can be calculated by performing a Lorentz boost of (jg, multiplying 
it by F(r, Pe) [see Eq.(2)] and integrating over Pg- Then, the scattering kernel is given by 



1 + 
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where D — 1 — iJ,V, Di — 1 — n'V, 7 = (1 — V^) and ^ = f2' • f2 is the cosine of scattering 
angle. In deriving the above equation we have chosen 
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(14) 



so that the distribution of electrons is normalized by a fixed electron density n as measured 
in the orthonormal frame associated with (1). 



3. THE METHOD OF SOLUTION 



3.1. Sepeiration variables 



As long as the ejected low energy photons satisfy z' — /ii/o/"^eC^7 <^ 1, the integration 
over incoming frequencies vq is trivially implemented provided that the explicit function of 
A^(r, fo, i^, O) is known. Thus, we need to describe the main properties of Green function 
N{r^ uo, u, ft) in a situation when the low-energy photons are injected into atmosphere with 
the bulk motion. 

The power-law part of the spectrum (Sunyaev & Titarchuk 1980, TL95) occurs at 
frequencies lower than that of Wien cut-off {E < E^, where Eg is the average electron 
energy). In this regime the energy change due to the recoil effect of the electron can be 
neglected in comparison with the Doppler shift of the photon. Hence we can drop the third 
term in parenthesis and the term ^ — 1 of the delta-function argument in the scattering 
kernel (13) transforming that into the classical Thomson scattering kernel (cf. TL95 and 
Gieseler & Kirk 1997). 

Now we seek the solution of the Boltzmann equation (12) with the aforementioned 
simplifications, in the form 



('+")j(r,At). 



(15) 



Then we can formally get from (12) that 



i~oJ 
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dr r J d/i 




(16) 



Here the phase function -R(^) 




(17) 
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where ^' is the cosine of scattering angle between photon incoming and outgoing directions 
in the electron rest frame. 

The reduced integro-differential equation is two dimensional and it can be treated and 
solved much easer than the original equation (12). The whole problem is reduced to the 
eigenvalue problem for equation (16). We can not claim that the kinetic equation allows a 
power-law solution (15) unless first a is found and J{r,ii) is specified. 

In order to derive equation for the determination of a spectral index we expand the 
phase function R{^) in series of Legendre polynomials (see also Sobolev 1975 and TL95) 



and 



R{0 = /^') + 2 E /^O COS - ^'), (18) 

m=l 

n 

p^{f^,fj.') = j2crpr{f^)pi^{f^') (19) 



C - C^^-^ (20) 
n + my. 



for m — 0, 1, 2, ... n. Since the phase function R{$,) is given by the series (Eq. 18) in 
cosm</7, the source function (the second term in brackets of the right hand side of Eq. 16), 
and J(r, O) can be expanded over cosm<p too. Under assumption of spherical symmetry for 
the source, we are interested in the zero-term of the expansion which satisfies the following 
equations 

£J°(r,/x) = -Kar + (a + 3)/x^] J°(r, /x) + (near)5°(r, /x), (21) 

where 



and the source function 

1 /•! 



B%T,^x) ^-J_y{fx,fi')Ar,fxW- (23) 



-14- 



There are two boundary conditions which our solution must satisfy. The first is that there 
is no scattered radiation outside of the atmosphere 

J°(0,//) = for fx<0. (24a) 

The second boundary condition is that we have an absorptive boundary at radius 

J^{r„fj,)^0 for //>0. (246) 

Thus the whole problem is reduced to the standard radiative transfer problem for the 
space part of the solution J(r, f2). Inversion of the differential operator £ of the left hand 
side of equation (21) leads to the integral equation for B^{r,iJ,) 

X "" exp{-T[nn,r'{r,i^'),i^']}B''{r',i^')dT 
Jo 

X exp{-T[r^,r'{r,i,'),i,']}B\r',i,')dT, (25) 

Jo 

where T{rbn, r, fi) is the optical path along characteristic curve of the differential operator £ 
determined by the initial point r, ii' toward the boundary radius r^jj (r^^ = r^, and oo for 
the inner, and outer boundaries, respectively). 

The phase function component p°(/x, /x') entered in equations (21) and (23) is 
determined by the sum 

p^{li,lx')^Y^C,P,{ii)PM. (26) 

i=0 

Thus, we can present the source function also as a sum: 

n „i 

B%r,i,) = Y.QP,{i,) P,{ii')J%r,ii')dii'. (27) 
i=o -^"^ 
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This form of the source function is used for the solution of the boundary problem (21-24) 
by the iteration method (e.g. Sunyaev & Titarchuk 1985). In order to proceed with the 
iteration method one has to assume some initial field distribution (in terms of the intensity 
J°) and then to calculate in accordance to Eq. (27) which is followed by the solution of 
the differential equation (21). 

This iteration formalism is identical to the integral-equation formalism in which 

S°(r,/.) = ^QP,(/.)SO(r) (28) 

i=0 

where the set of Bf{r), components of the source function B{r, fl), is determined by the 
system of the integral equations (compare with TL95) 

xj^ exp{-T[rbn,r'{r,i^'),i^']}B'j{r')dT 

+ f Piifi')Pj{fi')d(i' 
Jo 

xj^ eM-T[nn,r'(r,fx/),ii']}B^ir')dT]. (29) 

Thus the eigenvalue problem Eqs. (21-24), can be reduced to a eigenproblem for a system 
of integral equations (29) where the optical paths T{rbn, r, /i) and the expansion coefficients 
of the phase function Ci depend on the spectral index a as a parameter. In other words, 
one has to find the values of a which guarantee the existence of the nontrivial solution of 
equations (29). In §3.2 and Appendix B we shall proceed with the numerical solution of the 
eigenvalue problem by presenting the bulk motion phase function Rb{Cb) in the degenerated 
form (cf. Eq. 26). 

Now it is worth noting that in the case of the pure thermal motion in the isothermal 
plasma cloud the problem is substantially simplified. The source function B^{r,ii) can be 
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replaced by its zeroth moment -BQ(r) (TL95) which guarantees the accuracy of the spectral 
index determination better than 10% in the worst cases (Giesler & Kirk 1997). 

For example, the equation for the zeroth moment ^^(r) reads 

+ d^i' eM-T[rbn,r'ir,i^f),i^']}B^oir')dT]. (30) 

Jo Jo 

where Co is the zero-moment of the phase function. 



3.2. Photon trajectories and the characteristics of the space operator £ 

The characteristics of the differential operator i are determined by the following 
differential equation 

1 



dx ^ d[hi{i - (31) 



where x — r/rg is a dimensionless radius. The integral curves of this equation (the 
characteristic curves) are given by 

x(l-^^Y/' _ Xo(l - /i^)V^ _ ^ 

(l-a;-i)V2 (l-Xo-^)V2 ^' ^''^> 

where p is an impact parameter at infinity, p can also be determined at a given point in a 
characteristic by the cosine of an angle between the tangent to and the radius vector to the 
point and by the given point position xq. 

In the the flat geometry, the characteristics are just straight lines 

x{l-i,Y'=P, (33) 
where an impact parameter p is the distance of a given point to the center. 
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We can resolve equation (32) with respect of // to get 



/. = ±(l-pVy^) 



(34) 



where y — x^l"^ j {x — 1)^/^. The graph of y as a function of x is presented in Fig. 1 which 
allows to comprehend the possible range of radii for the given impact parameter p through 
the inequality p < y. For example, if p < V6.75, then the photon can escape from the 
inner boundary (the black hole horizon) toward the observer or vice versa all photons going 
toward the horizon having these impact parameters are gravitationally attracted by the 
black hole. However, if p > \/6.75, then the finite trajectories are possible with the radius 
range between 1 < x < 1.5, or the infinite trajectories with p < y{x) {x is always more than 



We are assuming a free fall for the background flow where the bulk velocity of the 



rh — M /Me, M is mass accretion rate and Me = Le/c^ — 'inGMnip/aTC is the Eddington 
accretion rate. 

For the cold converging inflow {kT^ — keV) the electron distribution is the 
delta- function F(r, Pg) = S{v — Vb) deflned in the velocity phase space in the way that 



1.5). 



3.3. Spectral index determination 



infalling plasma is given by v{r) = c{rs/rY^'^. In the kinetic equations (12, 16) the density 
n is measured in the local rest frame of the flow and it is n = m{r gj r)^^"^ / {^rax) ■ Here 




In this case the phase function is 




(35) 
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where subscript "b" is related with the bulk velocity direction (the case of arbitrary 
temperature will be considered elsewhere). In the case of zero temperature the directions 
of incoming and outcoming photons are related. Our goal is to find the nontrivial solution 
J^{r,fi) of this homogeneous problem and the appropriate spectral index a for which this 
solution exists. This problem can be solved by the iteration method which involves the 
integration of the differential equation (21) with the given boundary conditions (24) along 
the characteristics (32) by using Runge-Kutta's method. 

The integration starts from the internal or the outer boundary depending on the 
particular impact parameter p (Eq. 32). In turn which is determined by the dimensionless 
radius x {x = r /rg) and the cosine /i of the angle between the photon direction and the 
radius vector at the given point x. 

If /i is positive at the x then the photon trajectory (the characteristics) can start at the 
inner boundary [ if x < 1.5 and p < x^/^/(x — 1)^/^ or if a; > 1.5 and p < V6.75] or at the 
outer boundary (if x > 1.5 and p > V6.75). 

All cases can be understood from Fig. 1. The trajectories with the given p are related 
to the parallel lines to the X-axis, y = p- These lines start at x = 1 or at infinity. For 
example if they start at infinity (i.e. having negative /x) and p > -\/6.75 they must have the 
turning point with /i = 0. Thus they must pass through the point with radius where 
p = xj /{Xi, — 1). At this point the cosine changes sign from — to + and after that the 
trajectory enters through the point with radius x at the positive angle 9 = cos~^ /i. 

EDITOR: PLACE FIGURE |1| HERE. 

If the trajectory starts at x = 1 , (i.e. having positive /x) and p < V6.75 the parallel 
line y = p has no turning point. 
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If /i is negative at x, the trajectories starting at the internal boundary (having positive 
cosine /i) have to pass through the turning point /i = (changing the cosine sign) at radius 
Xi, where p = xV^ /{x^, — 1). However if the trajectory starts at the outer boundary has no 
the turning points the cosine /i is always negative along the trajectory. 

This space integration is followed by integration J°(r, yu') over the angular variable /i' in 
Eq. (23). As the initial distribution for B^{r,fi) or J^{r,fi) we can choose, for example, the 
uniform one. We use the Gaussian integration to calculate of Eq. (23) (see e.g. Abramovitz 
& Stegan 1970 for details of the methods). After quite a few iterations the iterative process 
converges and it produces the eigenfunction source distribution B^{r,fi). The number of 
iterations n is related to the average number of scatterings which the soft photons undergo 
to transform into the hard ones (e.g. Titarchuk 1994). It is determined by the Thomson 
optical depth of the bulk motion atmosphere = TT{rs). For the cold atmosphere (Tg = 0) 
the iteration number n ~ 2rf,. The convergence of the process can be done only with the 
proper choice of the value of a spectral index a. 



4. RESULTS OF CALCULATIONS AND DISCUSSION 

Fig. 2 presents the results of the calculations of the spectral indices as a function of 
mass accretion rates. It is clearly seen that the spectral index is a weak function of mass 
accretion rate in a wide range of m = 3 — 10. The asymptotic value of the spectral index 
for high mass accretion rate is 1.75 which is between a = 2 that was found by Blandford & 
Payne 1981 for the infinite medium and a ^ lA which was found by TMK96 for the finite 
bulk motion atmosphere. 

EDITOR: PLACE FIGURE | HERE. 
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The latter two results are obtained in the nonrelativistic Fokker-Planck approximation. 
We see that the efficiency of the hard photon production in the cold bulk motion atmosphere 
is quite low. This is not the case if the plasma temperature is of order of a few keV or 
higher. The coupling effect between the bulk and local Maxwellian motion occurs when 
the bulk motion velocity is very close to the speed of light, i.e. when the matter is very 
close to the horizon. The upscattering effect increases significantly in the latter case. In the 
regime of the relativistic bulk motion the electron distribution (2) has a sharp maximum 
at = TT and V = c. In the vicinity of the maximum the distribution is characterized by 
the exponential shape F{r.Pe) oc exp(— /57fe/27) (see also Eq. 5). More results and details 
regarding the relativistic coupling would be presented elsewhere. 

EDITOR: PLACE FIGURE | HERE. 

As an example, in Fig. 3 we demonstrate the zeroth moment of the source function 
distribution (the hard photon production). It is seen there that the distribution has a 
strong peak around 2 r^. This means that the vicinity of the black hole is a place where the 
hard photons are produced by upscattering of the soft photons off the converging electrons. 

Our calculations were made under the assumption of the free fall velocity proffie. Since 
the energy gain due to the bulk motion Comptonization is not bigger than factor 3 (if the 
spectral indices are higher than 1.5, TMK97), it follows that we can safely neglect the 
effects of the radiation force in our calculations if the injected photon flux in the converging 
inflow is of order of a few percent of the Eddington luminosity. 

The assumption of Thomson scattering accepted in our solution restricts the relevant 
energy range to E < mf,c^. Our approach cannot determine accurately the exact position of 
the high energy cutoff which is formed due the downscattering of the very energetic photons 
in the bulk motion electron atmosphere. The additional efforts are required to conflrm 
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the qualitative estimates of the high energy cutoff position as of order mgC^ (TMK97). 
Laurent & Titarchuk 1997 (in preparation) by using Monte Carlo calculations checked and 
confirmed our results for the spectral indices and the TMK97 estimates of the high energy 
cutoff position. Futhermore, they found the prominent spectral features at energies ^ 400 
keV. 

As a conclusion we would like to point out the definitive (according to our model) 
difference between black holes and neutron stars, as it can be ascertained in their spectral 
properties while in their soft states, when their luminosity is dominated by the quasithermal, 
soft, component: In the black hole case there should always be an additional steep power 
law high energy tail extending to energies ~ mgC^. This component should be absent in 
neutron star systems, because the effect of the bulk motion is suppressed by the radiation 
pressure in this case. 

We presented the full relativistic formalism and solved semi-analytically the kinetic 
equation by using TL95 eigenfunction method (see also Giesler & Kirk 1997) in the case 
of plasma inf ailing radially into a compact object with a soft source of input photons. We 
found that the converging flow has crucial effects on the emergent spectrum for moderately 

super- Eddington mass accretion rates. 

Our power law spectra can be applicable for the explanation of the observational 
situations in black hole candidate sources. 
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A. GENERAL RELATIVISTIC RADIATIVE KINETIC FORMALISM 

Let u be an arbitrary time like unit vector (u • u = u°'Ua = — 1, u° > 0) at some given 
space time point x, and let dix, d2X, d^x be the three arbitrary displacement vectors (with 
components dix", etc) which span an element of hypersurface orthogonal to u. By using 
orthogonality of two vectors ux and 

to the element of the hypersurface one can get the invariant volume element orthogonal to 
the unit vector u as follows: 

dV = ^/{^exapydix''d2X^d3x\ {A2) 

where g — det\gap\ , QajS is the metric tenzor (see Eq. 1) and e^a/j-y is the Levi-Civita 
alternating symbol eoi23 = +1- Since is independent of u we may orient dV to make 
vP = 1, = 0. Thus in a local Minkowski frame it would be 

dV = d^x. {A3) 

Similarly, at the point p (p • p = 0) the 3-surface element of the zero-mass shell spanned by 
three displacement vectors dip, d2P, d^p is presented by 

{dP)px = ^/{^)ea0,xdip''d2P^dsp\ {A4) 

Thus the invariant volume dP is 

; dip'd2P>d^p^ 

dP = \/{-g)eijk — ■ (^5) 

because all products of the right hand side of equation (A4) with p^ for A 7^ are equal 
to zero [this follows from equality p • p = 0, i.e. dp^ — {gijp' dp'^) / goaP'^^jf^i^ and the fact 
that a numbers of the even and odd transpositions of the combination a(3^ are equal). 
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When local Minskowian coordinates are used, with dip, d2P, d^p tangent to the respective 
coordinate lines, this reduces to the familiar Lorentz invariant 3— surface element in the 
momentum space, (Pp/E (e.g. Landau & Lifshitz 1971). If we introduce in the zero-shell 
spherical coordinates {E, 9, (p) defined (in a local Minkowski frame) through relationships 

p^ — E cos 9, p^ — E sin 9 cos ip, p^ — E sin 9 sin (f (A6) 

where £^ is a photon energy, then we obtain 

dP = pdEdn. {Al) 

In any coordinate system one can introduce the invariant infinitesimal volume element over 
the photon space, erected above the space time point under consideration. In other words 
one can introduce an orthonormal tetrad in each point of space-time and get the same 
result as equation (A7). 

The number of photon world lines crossing an infinitesimal "3-area" dTi\ (or dV) at x, 
with 4-momenta in the range dP is 

dN = N{x, p) = iV(x, p) (-P ■ vL)dPdV (A8) 

where the inclusion of the projector (— p • u) takes care of perpendicular counting of photon 
states along dUx or dV. 

The above somewhat abstract but fully covariant definition of the photon distribution 
function N incorporates, in fact, the familiar density interpretation of the classical 
Boltzmann distribution function defined over the six-dimension classical phase space of the 
photon gas. One simply has to take dV along an arbitrary t — const (see Eq. A3) spacelike 
hypersurface equipped with orthogonal coordinates, to recover from (A8) and (A7) that 

dN = N{x, p)d^xd^P, {A9) 
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i.e. the familiar phase space density interpretation of the distribution function. Comparing 
(A9) with the classical definition, 

dN = f{x, E, n)dE dQd^x, {AW) 

(Q is the unit 3-vector in the direction of the beam) one can obtain that 

N{x,p)^E-''f{x,E,n) (All) 

in the local proper frame. Since is a scalar, this gives its value in any other frame as 
well; it can be used to infer the transformation properties of / and the specific intensity 
I^^Ef = E^N. 

The general relativistic transfer equation essentially describes the rate of change of the 
above defined number of world lines dN as it is pushed along an arbitrary null spacetime 
geodesic. In the presence of external fields, the number of photon world lines centered 
around p, is altered. The alteration depends upon the specific interaction the photon 
field is subjected to. In accordance with the relativistic form of Liouville's Theorem (e.g. 
Lindquist 1966) (— p • u)dVdP remains invariant along the given set of trajectories. Hence 
the change in the number of world fines within (— p • u)dVdPdT — dW dP, where r is some 
parameter changing along geodesies (Landau & Lifshitz 1971) , is 



which is simply proportional to the change in N . Thus the transfer equation takes the 
following form 

where a product >S'(A^) dW dP denotes the change dN in all possible interactions taking 
place between the radiation field and external sources within dW . F^^ are the Christoffel 
symbols. 
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In the present paper we shall be interested exclusively in the Compton scattering of 
the radiation field of a locally Maxwellian distribution of electrons representing the plasma 
inflowing into black hole. To deal with the general covariant of the theory we shall perform 
all the calculations relative to the local orthonormal frame or tetrad {ea{x)} (a = 0, 1, 2, 3) 
tied to the metric (1) (e.g. Landau & Lifshitz 1971). Relative to such a local orthonormal 
frame, an arbitrary photon momentum p may be represented as 

P = pX- (^14) 
If is the usual contravariant components of p in some coordinate system {x""} then 

P = p"e«. {A15) 
where {6^} is the induced coordinate basis. If we deflne 

ea = e^e^, and = {A16) 

then it is evident that 

p'' = e>" and = {A17) 

correspondingly. One can introduce the transformation which leaves the coordinate system 
unaltered, but express vectors in terms of their tetrad components, namely 

x"" = x". 

By using this transformation the radiative transfer equation can be rewritten in the form, 
which is very similar to Eq. (A13) 

dN dN u 



with 
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being the Ricci rotation coefficients (e.g. Landau & Lifshitz 1971). 



We introduce in the tangent space, at each event {t, r, 9, <^), the orthormal basis 
Bq = e^^^et, ei = e~^^^er 62 = r'^e^ 63 = (r sin 6')"^e^, 



{A21) 



where e^/^ = / [see eq. (1)]. Relative to this local orthormal frame an arbitrary photon 
momentum can be represented as 



p= = (£',£' cos ^, £" sin 6* cos sin 6* sin 99) 



{A22) 



where E is the physical photon energy as measured by an orthonormal observer, and 9, (p 
are local polar-like phase space coordinates. The angle 9 is measured relative to the outward 
radial direction at the point under consideration. For our problem we shall assume the 
steady state conditions for the accretion and radiation field, respectively, implying that N is 
time-independent and spherically symmetric. Furthermore, A'^ should be axially symmetric 
around the radial direction, which results in = N{r, E, fi) where n = cos9. 

Thus the coefficients e° and are given by 



{A23) 



and 



et^S^r, e2^52r sin 9. 
Here Sj is Kronecker delta (= if i 7^ j; = 1 if i = j). 



(A24) 



By using these formulas one can calculate the rotational coefficients 7^^ [see eq. (A20)] 
and get the relativistic radiative transfer equation in the form 



-A/2 



dN (1 - fi"^) dN 



+ e 



-A/2 



A' 



ON , 



E ' 



(^25) 
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The terms in the first bracket are obtained by using formulas (A24) for e|j^°) and e^^^) 
the fact that the derivative 

dO 

The first term in the second bracket is related with the Ricci coefficient 7qq and the second 
term there appears as a result of calculation of product ■Jooip^y {ON / dp^) . We also use the 
fact that the partial derivative 

ON l-n'^dN ON 
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B. NUMERICAL CALCULATION OF THE SOURCE FUNCTION 

We use the iteration method to solve the radiative transfer problem Eqs. (21-24) which 
involves the calculation of the source function B^. In order to calculate the source function 
we transform it into the degenerate form by using the phase function formula (35). We have 

5°(r,/.) = ^ j\\ii,ti')j\r,ti')dti' = {h + l2 + h + h) + (/3 + /s)/^ + {h - 17)1^'. (Bl) 

The following integrals Jj are computed by using the Gaussian integration formula (see for 
details of the Gaussian integration method, e.g. Stegan & Abramovitz 1966). 



^ /-I ^^^^ 

^ lb -^-1 

^ ^^^^ 

^ lb ^b 
° lb ^b 
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Fig. 1. — Plot of the photon trajectory phase space. The phase space function y versus the 
dimensionless radius x = r/r^. 

Fig. 2. — Plot of the energy spectral index (photo index minus one) versus the total mass 
accretion rate. 

Fig. 3. — Plot of the source function distribution, Eq. (Bl) in arbitrary units versus the 
dimensionless radius x = r/r^ for the dimensionless mass accretion rate m = 4 and = 0. 
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